A three‐dimensional Cartesian cut cell method for incompressible viscous flow with irregular domains

A three‐dimensional Cartesion cut cell method is presented for the simulations of incompressible viscous flows with irregular domains. A new model (referred to as ‘6+N’ model) is proposed to describe arbitrarily shaped cut cells and treat all the cells as polyhedrons with 6+N faces. The finite volume discretization of the Navier–Stokes equation is then implemented by using the ‘6+N’ model to separate the surface flux integrals into two parts, that is, the fluxes through the basic face of the hexahedron and those through the cutting surfaces. The previously proposed Kitta Cube algorithm and volume computer‐aided design platform (J. Comput. Aided. Des. 2005; 37(4): 1509–1520. Doi:10.1016/j.cad.2005.03.006) are adopted to generate cut cells and provide shape data and physical attributes for the numerical analysis. A modified SIMPLE‐based smoothing pressure correction scheme is applied to suppress checkerboard pressure oscillations caused by the collocated arrangement of velocities and pressure. The calculation accuracy of the numerical method expressed by L₁ and L _∞ norm errors is first demonstrated by the simulation of a pipe flow. Then its feasibility, efficiency, and potential in engineering applications are verified by applying it to solve natural convections between concentric spheres and between eccentric spheres. The heat transfer patterns in eccentric spheres are also obtained by using the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.     

求解HJB方程的拟变分不等式组的迭代法

本文对HJB方程的拟变分不等式组提出一种迭代算法,并给出此算法在一定的条件下的单调性定理和证明,数值试验表明此法有效的.     

Implementation of ADER scheme for a bore on an unsaturated permeable slope

The paper details the implementation of the Godunov‐type finite volume Arbitrary high order schemes using Derivatives (ADER) scheme for the case of a large source term in the continuity equation of the nonlinear shallow water equations. The particular application is the movement of a bore on a highly permeable slope. The large source term is caused by the infiltration into the initially unsaturated slope material. Infiltration is modelled as vertical downwards piston‐like flow with Forchheimer quadratic parameterisation of the resistance law. The corresponding ODE is solved using the fourth‐order Runge–Kutta method. The surface and subsurface flow models have been tested by comparison with analytical solutions. Example predictions of surface bore propagation and wetting front propagation are presented for a range of slope permeabilities. The effects of permeability on bore run‐up, water depths and velocities are illustrated. The ADER scheme is capable of handling the source term, including the extreme case when this term dominates the volume balance. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation and convergence analysis for a system of nonlinear singularly perturbed differential equations arising in population dynamics

In this article, we consider a system of nonlinear singularly perturbed differential equations with two different parameters. To solve this system, we develop a weighted monotone hybrid scheme on a nonuniform mesh. The proposed scheme is a combination of the midpoint scheme and the upwind scheme involving the weight parameters. The weight parameters enable the method to switch automatically from the midpoint scheme to the upwind scheme as the nodal points start moving from the inner region to the outer region. The nonuniform mesh in particular the adaptive grid is constructed using the idea of equidistributing a positive monitor function involving the solution gradient. The method is shown to be second order convergent with respect to the small parameters. Numerical experiments are presented to show the robustness of the proposed scheme and indicate that the estimate is optimal.     

A new coding scheme of QC-LDPC codes for optical transmission systems

Based on Galois Field (GF(q)) multiplicative group, a new coding scheme for Quasi-Cyclic Low-Density Parity-Check (QC-LDPC) codes is proposed, and the new coding scheme has some advantages such as the simpler construction, the easier implementation encoding, the lower complexity of the encoding and decoding, the more flexible adjustment of the code length as well as the code rate and so forth. Under the condition of considering the characteristics of optical transmission systems, an irregular QC-LDPC (3843,3603) code to be suitable for optical transmission systems is constructed by applying the proposed new coding scheme. The simulation result shows that the net coding gain (NCG) of the irregular QC-LDPC (3843,3603) code is respectively improved 2.14 dB, 1.19 dB, 0.24 dB and 0.14 dB more than those of the classic RS (255,239) code in ITU-T G.975, the LDPC (32640,30592) code in ITU-T G.975.1, the regular SCG-LDPC (3969,3720) code constructed by the Systematically Constructed Gallager (SCG) coding scheme and the regular QC-LDPC (4221,3956) code at the bit error rate (BER) of 10-8. Furthermore, all the five codes have the same code rate of 93.7%. Therefore, the irregular QC-LDPC (3843,3603) code constructed by the proposed new coding scheme has the more excellent error-correction performance and can be better suitable for optical transmission systems.     

High-resolution semi-discrete Hermite central-upwind scheme for multidimensional Hamilton–Jacobi equations

We introduce a high resolution fifth-order semi-discrete Hermite central-upwind scheme for multidimensional Hamilton–Jacobi equations. The numerical fluxes of the scheme are constructed by Hermite polynomials which can be obtained by using the short-time assignment of the first derivatives. The extensions of the proposed semi-discrete Hermite central-upwind scheme to multidimensional cases are straightforward. The accuracy, efficiency and stability properties of our schemes are finally demonstrated via a variety of numerical examples.     

An anti‐diffusive volume of fluid method for interfacial fluid flows

A volume of fluid (VOF) method is developed combining a first‐order limited downwind scheme with higher order accurate schemes. The method is characterized by retaining a sharp fluid interface and a reduction in numerical diffusion near the interface, but avoids complicated geometrical reconstruction as occurs in most volume tracing algorithms. To demonstrate the accuracy and robustness of the method, a selection of numerical experiments are presented involving a pure advection problem, a water wave impact caused by a dam breaking and liquid sloshing in a partially filled tank. Copyright © 2011 John Wiley & Sons, Ltd.     

A 1D Roe decomposition of a system of equations governing wave height transformation

A Roe‐type decomposition for a system of equations governing onshore/offshore wave transformation in coastal waters is derived. The equation set approximated pertains to coastal waters prior to wave breaking, and is based on depth‐averaging and time‐averaging of the Euler equations. The equations are those used in many commercial codes for simulation of wave height and wave‐averaged currents. This novel approach uses a combination of some standard Roe averages, together with physical reasoning and power series expansions to derive a Roe‐averaged Jacobian (with real, linearly independent eigenvectors) and ensures conservation, and thereby effects the decomposition. It is shown that the resulting derived Roe‐averaged quantities are accurate to a high degree, by comparing them with their analytical equivalents for a wide range of nondimensional water depths and slopes likely to be encountered in coastal problems. Numerical tests of time‐invariant wave height transformation and wave group propagation are undertaken; these indicate good performance of the scheme in practice. Copyright © 2011 John Wiley & Sons, Ltd.